Non-Markovian dynamics of double quantum dot charge qubits due to acoustic phonons 
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We investigate the dynamics of a double quantum dot charge qubit which is coupled to piezoelectric acoustic 
phonons, appropriate for GaAs heterostructures. At low temperatures, the phonon bath induces a non-Markovian 
dynamical behavior of the oscillations between the two charge states of the double quantum dot. Upon applying 
the numerically exact quasiadiabatic propagator path-integral scheme, the reduced density matrix of the charge 
qubit is calculated, thereby avoiding the Born-Markov approximation. This allows a systematic study of the 
dependence of the Q-factor on the lattice temperature, on the size of the quantum dots, as well as on the interdot 
coupling. We calculate the Q-factor for a recently realized experimental setup and find that it is two orders 
of magnitudes larger than the measured value, indicating that the decoherence due to phonons is a subordinate 
mechanism. 

PACS numbers: 03.67.Lx, 03.65.Yz, 73.63.Kv, 63.20.Kr 



I. INTRODUCTION 

Various candidates for realizing building blocks of quan- 
tum information processors with nanoscale solid state struc- 
tures have been proposed and also partially realized in ground- 
breaking experiments. An important class of proposals con- 
sists of using the charge degree of freedom in semiconducting 
double quantum dots (DQDs)i£ to realize a quantum mechan- 
ical two-state system or quantum bit (qubit) i^ii^ Thereby, the 
logical states |0) and |1) are given by the low-energy charge 
states of the DQD with one excess electron either on the left 
or the right dot. The transition between these states occurs via 
tunneling between the two dots. A significant advantage of 
the charge qubit is its direct controllability via external volt- 
age sources. In recent experiments, the coherent manipulation 
of charge states of DQDs has been achieved. 7,8,9 Thereby, the 
control over the bias and the gate voltages permits to reliably 
tune the DQD in the Coulomb blockade regime to the required 
transition line (1,0) <-> (0, 1). The large charging energies 
suppress leakage to energetically higher lying many-electron 
states. 

A central impediment for coherent quantum information 
processes is decoherence and dissipation, see Ref. [T(J for a 
particular example. Since charge qubits are rather easily ad- 
dressable from outside, they are, in turn, also rather fragile 
for various decoherence mechanisms from the environment. 
In order to achieve a thorough understanding of the role of 
the various decoherence mechanisms, one has to rely on re- 
alistic and precise model calculations which allow to sort out 
the different contributions. Among them, decoherence due 
to acoustic phonons is one possible candidate. This mecha- 
nism has been studied previously upon using approaches rely- 
ing on the Born-Markov approximation^iiiii^ The decoher- 
ence rates obtained in these studies were one or two orders of 
magnitude smaller than those determined experimentally for 
GaAs DQDs. 7,8 However, no exact treatment was available to 
validate the Born-Markov approximation in the DQD context. 
This is particularly interesting in view of realistic DQD de- 
signs since geometrical form factors tailor a rather peculiar 



spectral density J(u>) of the phonon environment. For piezo- 
electric phonons, it has a high-frequency tail decaying only 
algebraically with!! J(w) oc uj^ 1 while it is superohmic at 
low frequencies with J(u>) oc ui 3 , with a crossover at frequen- 
cies in the regime of the tunneling amplitude of the DQD. As 
it has also been pointed out in Ref. 11, the use of the Born- 
Markov approximation is appropriate only at small tunneling 
amplitudes. However, is is expected to become increasingly 
unreliable at DQD with larger interdot tunneling amplitudes. 
The motivation to find exact reference solutions for this prob- 
lems also stems from the fact that the Born-Markov approx- 
imation has led to the conclusion that other mechanisms of 
decoherence, as for instance, background charge fluctuations 
and electromagnetic noise from voltage fluctuations, would be 
the dominating source of decoherence while phonon decoher- 
ence should be a negligible contribution. In view of further 
design optimization for real devices, exact results are clearly 
desirable. 

In this paper, we study phonon decoherence using the 
same model as in Ref. lll| as a basis for the numerically ex- 
act iterative quasiadiabatic propagator path integral (QUAPI) 
scheme piiiiii^ In particular, we avoid the Born-Markov ap- 
proximation which turns out to be problematic at larger tun- 
neling amplitudes A of the DQD. The decaying oscillations 
of DQD charge states allows us to extract the quality (or Q- 
) factor. We find that Q decreases with increasing tunneling 
amplitudes A for small A. In fact, the solution in this regime 
is accurately described by the weak-coupling result for a su- 
perohmic environment. 16 At large A, the Q-factor increases 
with increasing A and the non-Markovian corrections become 
noticeable changing even the order of magnitude of Q. In 
between these two regimes, a crossover occurs which is re- 
lated to the geometry of the DQD, see below. In particular, 
our numerical precise method allows to obtain results for the 
DQD device recently realized in experiments. 7,8 We find that 
the calculated value lies approximately two orders of magni- 
tude above the experimentally measured value. Indeed, this 
suggests that the phonon decoherence mechanism is subor- 
dinate and other mechanisms dominate decoherence in DQD 
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charge states which have been realized until now. In compar- 
ison to the results from the Born-Markov approximation, we 
find that the exact Q-factor is in general larger, indicating that 
the former overestimates phonon decoherence. 



II. MODEL FOR CHARGE QUBIT AND PHONON BATH 



with g q = AqF(q)(l — e _tq ' deH ). Notice that the phonons 
can propagate in all three dimensions and that the electron- 
phonon-coupling is not isotropic. 11 Having specified the cou- 
pling, we can now introduce the spectral density of the phonon 
bath which reads 



J M = ^2\9q\ 2 S{uJ ~ W q ) • 



(7) 



In order to solely investigate the role of the phonons in the 
decoherence processes, we start with the simplifying assump- 
tion that the DQD is isolated from the leads. The total Hamil- 
tonian is given by 6 - 16 



H = H s + Hb + Hsb- 



(1) 



Here, Hs is the Hamiltonian of the DQD which is modeled in 
the basis of the two localized charge states |0) = \L) = (1, 0) 
and |1) = \R) = (0, 1) as a symmetric two-level-system in 
pseudo-spin notation (cr, are the Pauli matrices) as 



(2) 



In this simple model, we have assumed that the external gate 
voltage is tuned such that the two charge states are energeti- 
cally nearly degenerate and that the excess electron can tun- 
nel between the two dots with the tunneling amplitude A. A 
possible energy bias between the two states could easily be 
included in the model but is not considered in this work. 
The Hamiltonian for the phonon bath is, as usual, given by 



Hi 



(3) 



with the phonon dispersion relation w q . 

The interaction part of the Hamiltonian is determined by 
the electron-phonon interaction and reads^AAiii 



H. 
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<^)(&q + 6-q). 



(4) 



Here, N% is the number of the excess electrons on the left/right 
(£ =L/R) dot and a q = A q e- iqR «F£(q) where R L = and 
Rr = de y are the position vectors of the two dots. As in Ref. 
Mil we have assumed that the center of the left dot is located 
at the origin of the coordinate system while the center of the 
right dot is placed on the y-axis with distance d, see Fig.^for 
a sketch. The two dots are assumed to have an equal radius of 
a and to be confined to the xy-plane in the underlying 2DEG. 
The coupling constants A q which depend on material param- 
eters and on the wave vector q will be specified below. To 
take into account a realistic dot geometry, we associate a form 
factor 
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to each dot, where (r) is the charge density distribution of 
the dot. In the following, the two dots are assumed to have 
identical form factors. With this the coupling Hamiltonian 
simplifies toiA*ii 



H S b = ^ ^ 5q (6 q + 5 -q) 



(6) 



As in Ref. [HI we specialize to linear acoustic phonons with 
velocity s (s w 5 x 10 3 m/s for GaAs) and linear dispersion 
w q = s|q|. Moreover, we only consider coupling to longitudi- 
nal piezoelectric phonons and neglect the contribution of the 
deformation potential which is justifiable at temperatures be- 
low 10 K for bulk GaAs material. 11,18 This yields the coupling 
constant |A q | 2 = <? p h7r 2 s 2 /(V|q|), where g p h is the dimen- 
sionless piezoelectric constant (g p h m 0.05 for GaAs) and V 
being the volume of the unit cell. Assuming that the charge 
density distribution on the dot is Gaussian in the icy-plane and 
localized in the z-direction, one finds for the spectral density 
of the bathii 

/•f/2 

J(w) = g ph io / d9 sin 0e- uj2a2sin2 6 /s2 



(8) 



The spectral density is shown in the inset of Fig.|2]for the case 
of GaAs and a dot radius of a = 60 nm. The low-frequency 
behavior is superohmic, i.e., J(uj — ► 0) « asow^w 3 with 
a so = 5ph^ 2 /(6a 2 ), while the high-frequency tail falls off 
algebraically like J(u> — > oo) oc l/io. The crossover be- 
tween these two limiting regimes occurs on a frequency scale 
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lj c = sj 'a = t~ . For GaAs and a typical dot radius of a = 60 
nm, we obtain oj c — 83 GHz corresponding to an energy of 
55 //eV. As we will see below, the typical tunneling ampli- 
tudes A are comparable to this energy scale. This indicates 
that the Markov approximation could become problematic in 
this transition region. The cut-off frequency uj c corresponds 
to the inverse of the time scale r c of the bath autocorrelation 
function which reads 16 



L(t) = L R (t) 
1 



HL/(t) 
dujJ{uj) 



coth cos tot 



i sin Lot 



(9) 



The autocorrelation function Lit) is shown in Fig.|2]for T = 
15 mK for a = 60 nm and d = 180 nm>iiiiS The algebraic 
decay at high frequencies determines the short time behavior 
of L{t); notice the finite slope of L^it) at t = 0, in contrast 
to zero slope of the usual exponential and Drude-shaped cut- 
off functions. 1 - The superohmic low frequency behavior of 
J(w) yields a fast decay of Lit) at long times. The often 
used Born-Markov approximation corresponds to replacing 
the strongly peaked real part Ln{t) by a ^-function with the 
corresponding weight while the imaginary part Liit) is often 
neglected. Moreover, we note that the dimensionless damp- 
ing constant given by the prefactor ago of the low-frequency 
superohmic limit assumes the value ago = 0.075 for the pa- 
rameters quoted above. Thus, our results at small A can be 
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compared to those from a weak-coupling approach in terms 
of real-time path-integrals for the spin-boson model with a 
superohmic environmentpA& 

IIL THE QUASIADIABATIC PATH-INTEGRAL 
PROPAGATOR (QUAPI) 

The dynamics of the qubit is characterized by the time evo- 
lution of the reduced density matrix p(t) which is obtained 
after tracing out the bath degrees of freedom, i.e., 

P (t) = tT B U{t,0)W(0)U- 1 {t,0), 

U(t,0) = Tcxpj-^ dt'H^ . (10) 

Here, U(t, 0) denotes the propagator of the full system plus 
bath and T denotes the time-ordering operator. W(0) is the 
full density operator at initial time set at t = 0. We as- 
sume standard factorizing initial conditions, 1 * i.e., W(0) oc 
p(0) exp[—Hs I (fcsT)], where the bath is at thermal equilib- 
rium at temperature T and the system is prepared according 
to p(0). Throughout this work, we always start from the qubit 
state p(0) = \L){L\. 

The technique which we use to calculate p(t) is the well es- 
tablished iterative tensor multiplication algorithm derived for 
the quasiadiabatic propagator path integral (QUAPI),— It is 
a numerically exact algorithm, as also, for instance, the real- 
time quantum Monte Carlo method is. 20 It has been success- 
fully tested and adopted in various problems of open quantum 
systems, with and without external driving>12iliil4 For details 
of this algorithm, we refer to previous worki2J2i£ and here 
only address the ingredients which are important to our calcu- 
lations. 

The algorithm relies on the symmetric Trotter splitting of 
the short-time propagator U(tk+i,tk) into system (Hs) and 
bath (Hb) parts on a time slice of length At. The bath dynam- 
ics can be solved exactly yielding a Feynman- Vernon-type in- 
fluence functional 16 while the system is propagated exactly 
by solving the Schrodinger equation. The neglect of higher 
order terms of the full propagator causes an error of the order 
of At 3 . The Trotter splitting allows to calculate an approx- 
imate value of the true result for the observable of interest, 
with an error depending on At. As shown in Ref. EI this 
error vanishes proportional to At 2 as At — > 0. Thus, by ex- 
trapolation of the results for At — > 0, the Trotter error can 
be eliminated completely (the details of this procedure will be 
discussed elsewhere 2 ^). 

The interaction of the system with the bath induces inter- 
nal transitions in both and creates intercorrelations between 
them (memory); the latter are described by the autocorrelation 
function L(t) given in Eq. For the phonon bath, these cor- 
relations decay on a time scale given by the correlation time 
r c . This motivates us to neglect such long-time correlations 
beyond a memory time r mom and to break up the influence 
kernels into pieces of total length r mcm = K At, where K 
denotes the number of time steps over which the memory is 
fully taken into account. 



The two strategies mentioned above are countercurrent: For 
the Trotter splitting, a small time step At is desirable, thus de- 
creasing r mom for a fixed K. On the other hand, a large r mcm 
is wanted in order to take a long memory range into account. 
Nevertheless, converged results can be obtained in the follow- 
ing way: (i) We choose T mem such that we include all the rele- 
vant parts of the correlation function L(t). Quantitatively, this 
can be done by requiring L i?(r mcm ) / L r(0) < 10~ 3 in the 
asymptotic regime, (ii) We choose K such that the resulting 
At = r mcm /K is small enough to ensure that we are in the 
regime which allows extrapolation for At — > (see above). 
This procedure allows to completely eliminate the Trotter er- 
ror. We note that the choice of K is limited to a maximum of 
K = 12 for reasonable numerical practicability on a personal 
computer. 



IV. THE Q-FACTOR OF COHERENT CHARGE 
OSCILLATIONS 

We have adjusted the QUAPI algorithm to the phonon spec- 
tral density and determine the charge population pn{t) 
of the left dot as a function of time for the initial condition 
Pn(O) = 1. A typical example of the resulting damped oscil- 
latory behavior is shown in Fig.[3]where the exact QUAPI re- 
sults are depicted by the symbols (♦) for the case A = 27peW, 
T = 15 mK and K = 12. To extract the damping rate 7 and 
the oscillation frequency ft, we fit an exponentially decaying 
cosine to the numerical data. The result of the fit is also shown 
in Fig.|5]as solid line. The fit yields SIquapi = 0.98cj c for the 
oscillation frequency and 7quapi = 5.27 x 10~ 3 w c for the 
decay constant. 

The ratio of 7 and ft fixes the quality factor according to 
the convention Q = Q/(ttj). Evaluating Q allows us to in- 
vestigate the dependence of the coherence of the charge oscil- 
lations on various experimental parameters. Figure [4] shows 
the Q-factor as a function of the tunneling amplitude A ob- 
tained with the numerically exact QUAPI algorithm (+). For 
illustration, we have used parameters from the experimental 
setup of Ref. d Since our model contains realistic assump- 
tions about the geometry of the DQD and the materials char- 
acteristics, the predictions of the QUAPI calculations in this 
case can be considered quite accurate. Two regimes exist: 
(i) At small A, Q decreases with increasing A (see below) 
while (ii) at large A, Q increases again with increasing A. 
The crossover in between occurs when the energy splitting 
between the qubit states coincides with the maximum of the 
spectral density J(uj), i.e., when 2A = ui c . Then, the damp- 
ing is most efficient and thus the rate 7 is maximal. Notice 
that realistic values for the tunneling amplitudes A fall in the 
range of the bath correlation frequency scale uj c (see upper 
axis of Fig. @J where A is scaled with respect to ui c ). 

Since the coupling to the phonons is rather weak (cf. the 
value of aso = 0.075 for the superohmic coupling constant at 
low frequencies), it is tempting to compare our exact results to 
simple, approximate analytical results obtained from real-time 
path-integral formalism in the regime of weak-coupling. 16 
However, we remark that one needs to have in mind that the 



4 



algebraic decay of the spectral density at large frequencies is 
by itself not taken into account while deriving the approximate 
result. Thus, only the superohmic low-frequency part can be 
expected to yield a reasonable agreement. The weak-coupling 
results are given byi& 

= 2A[1 - 2Reu(2iA)] 1/2 (11) 

and 

7 =Ij(2A)coth-^, (12) 
4 k B T 

where 

, , 1 f°° J(u) ( Hu \ 

^=u ^rw- 1 )' (l3) 

With this, we can compute the time evolution of charge pop- 
ulation of the left dot in the weak-coupling approximation 
which is shown in Fig. [5] (dashed line) for the same param- 
eters as above. The deviation from the exact result for the 
chosen parameters due to non-Markovian effects is striking 
and also illustrated by the numerical values £1 = 3.0ut c and 
7 = 5.6 x Kr 3 w c . 

The resulting Q-factor is shown in Fig. |4] (solid line, 
"wca"). Since the results in Eqs. (II It and (11 2t are obtained 
by linearization with respect to the damping constant g p ^, they 
are equivalent to a Born-Markovian approximation. We find 
noticeable deviations for intermediate and large values of A. 
From Eq. il Q it follows that the oscillation frequency fi is 
renormalized by the term involving u(2iA) stemming from 
one-phonon interblip correlations in the self-energy. 16 We find 
that at small A these corrections are negligible, while they in- 
creasingly become important at larger A, pointing to a non- 
Markovian dynamics in this regime. Note that the qualitative 
behavior of Q versus A is similar to that found in Ref. [TTJ 
within the Born-Markovian master equation approach. How- 
ever, the absolute numbers disagree. 

To understand the results at small tunneling amplitudes, we 
zoom into the small-A region and show the QUAPI results in 
Fig.|5](+, right ordinate scale) together with the results from 
the weak-coupling approximation (solid line) obtained from 
Eqs. il l\ and H2\ . Therefore, we have used the pure super- 
ohmic spectral density with an exponential cutoff with fre- 
quency il c 3> A, i.e, J(io) — asoui~ 2 LU 3 exp(— cj/f2 c ) with 
a so = 0.075 and fl c — 10cj c . We find a good agreement be- 
tween both results in the limit of small A, thus justifying that 
the dynamics at small A is Markovian. Nevertheless, note 
that even in this frequency range deviations between the exact 
and the approximate solution do appear. As a further check, 
we also compare the results of both the numerically exact and 
the analytical approach for Q for a pure Ohmic environmental 
spectral density J(lo) = auicxp(—ui/uj c ) for a = 0.05. Both 
the QUAPI (T) and the approximative weak-coupling solution 
(dashed line) coincide. The Q-factor decreases monotonically 
for decreasing A since the Ohmic low-frequency modes are 
not suppressed strongly enough in comparison with the super- 
ohmic case. Although the dynamics is Markovian for small 



A, non-Markovian corrections appear to be relevant at inter- 
mediate values of A. Because realistic values for the tunnel- 
ing splitting A are also of the order of u c , the role of the large 
frequency tail oc 1/uj still is not negligible and cannot be ap- 
propriately taken into account within a Born-Markovian mas- 
ter equation. 

The sensitivity of the Q-factor on the lattice temperature 
T of the GaAs host is shown in Fig.|6]for a small and a large 
tunneling amplitude A. Note that for the latter the correspond- 
ing Born-Markov results have been shown not to be reliable 
due to noticeable non-Markovian corrections which can be ob- 
served in our data. Another parameter which is adjustable in 
the DQD design is the dot radius a. The dependence of Q on 
a is depicted in Fig. for a large and a small tunneling am- 
plitude for a fixed ratio d/a = 3. We find that Q decreases 
monotonically with increasing a for small value of A, while it 
increases with increasing a for large value of A, qualitatively 
similar to the findings in Ref. 1 1 . Note, however, the size- 
able differences in the absolute numbers which are due to the 
non-Markovian corrections. 

Finally, we address the recent experiment of Hayashi et al.. 7 
These authors implemented a GaAs DQD in the bias-pulsing 
mode. During the pulse, the tunneling amplitude is constant. 
The bias pulse is ramped up on a time scale of 100 ps. The 
tunneling amplitude has been estimated as A = 5 /ieV. The 
lattice temperature was T = 20 mK, while the effective dot 
radius and the interdot distance have been estimated to be 
a « 50 nm and d ~ 225 nm, respectively. For this combina- 
tion of parameters, the QUAPI method yields Q = 352 which 
has to be compared with the experimental value of Q = 3. 
On the other hand, the weak-coupling solution of Eqs. (^J 
and H2i yields Q — 539. Two important conclusions can 
be drawn from these numbers. First, the Born-Markov ap- 
proximation underestimates decoherence in these systems and 
non-Markovian corrections are quite noticeable. Second, our 
exact result differs from the experimental value by a factor 
of roughly 100, indicating that the role of the phonons for 
the decoherence is indeed of minor importance to understand 
the reported experimental results. Other sources of decoher- 
ence like voltage fluctuations from the electromagnetic envi- 
ronment or background charge fluctuations exist and have to 
be taken into account for an accurate description for so far 
realized charge-based DQD qubits. 

However, our result go beyond this point: Our predicted 
exact values of the Q-factors represent a fundamental upper 
limit on the possibility to improve coherence of DQD charge 
qubits since phonon decoherence represents an intrinsic deco- 
herence mechanism which can hardly be overcome. This has 
to be noticed also in view of the DiVincenzo criteria 23 for a 
possible realization of a quantum computer. 



V. CONCLUSIONS 

To summarize, we have obtained numerically exact results 
for the Q-factor of the decaying charge oscillations in a double 
quantum dot upon using the real-time quasiadiabatic propaga- 
tor path integral (QUAPI). Realistic assumptions on the form 
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of the environmental phonon spectral density entered in our 
model via geometrical form factors and materials character- 
istics. No fitting parameters of any sort were utilized. We 
have investigated the quality {Q-) factor as a function of the 
tunneling splitting and have compared our results with those 
obtained from a weak-coupling approximation within an ana- 
lytical approach in terms of real-time path-integrals. We find 
that the regime of small tunneling amplitudes is appropriately 
covered by the Markovian description. However, at larger 
(but still realistic) values of the tunneling amplitude, non- 
Markovian corrections appears and are relevant. Moreover, 
we have determined the temperature dependence of Q-factor 
as well as its dependence on the dot radius. From a compar- 
ison with the result obtained in an experimental realization 
of a double quantum dot charge qubit, we find that the the- 
ory predicts a Q-factor two orders of magnitudes larger than 
the measured value. This leads to the conclusion that phonon 
decoherence is a subordinate mechanism in GaAs quantum 
dots realized with present state-of-the-art technology. Clearly, 
other forms of coupling to the external environment prevail 
for those already realized devices. However, our results also 
represent a fundamental upper limit to the coherence of DQD 



charge qubits which can hardly be overcome due to its intrin- 
sic nature. 

It would be interesting to extend the QUAPI calcula- 
tion to Si-based charge qubits, which have recently been 
implemented. 9 Since Si has inversion symmetry, piezoelectric 
phonons do not occur and one needs to include deformation 
potential phonons in the calculations. The decoherence in- 
duced by the electron-phonon coupling is likely to be smaller 
than in GaAs. However, for the deformation potential case, 
~ Thus, if on one hand we expect a strong super- 

ohmic behavior at low frequencies, on the other hand, a very 
slow decay of the spectral function at high frequencies should 
occur. That may provide an even stronger non-Markovian dy- 
namics and have implications for other implementations, such 
as Si:P donor-based charge qubits. 24 - 25 - 26 
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FIG. 1: Sketch of the geometry of the double quantum dot charge 
qubit and the various angles of the phonon propagation. 
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FIG. 3: Time evolution of the charge population of the left dot for 
the case A = 27/ieV and T = 15 mK. Shown are the exact results 
obtained from QUAPI (♦) and the result of a fit of an exponentially 
decaying cosine to pn(i) (solid line). The dashed line indicates the 
result of a weak-coupling (Born-Markov) approximation using Eqs. 
1 1 1 i and 1121 for the same parameters. Deviations from the exact 
results are striking. The remaining parameters are the same as in Fig. 
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FIG. 2: The bath autocorrelation function (response function) 
L(t) = Ln(t) + iLi(t) for the spectral density J(w) (inset) of the 
phonon bath for the case of GaAs with (? p h = 0.05, s — 5 x 10 3 
m/s, a dot radius of a — 60 nm, and interdot distance d = 180 nm. 
Temperature is T — 15 mK. The dashed lines in the inset mark the 



superohmic limit J(uj — > 0) = asoi^c w 
the high-frequency limit J(cu — ► oo) oc 1/u 



at low frequencies and 
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FIG. 4: The quality factor Q of the coherent charge oscillations as a 
function of the tunneling amplitude A in natural units (lower scale) 
and scaled with respect to lo c (upper scale). The data shown are 
obtained with QUAPI (+) using the phonon spectral density (Eq. ^8) 
and with the weak-coupling approximation (wca) of Eqs. i 1 H and 
1121 . The parameters are the same as in Fig.|2| 
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FIG. 5: Zoom into the small-A region. Shown are the QUAPI re- 
sults for the phonon (+) and the pure Ohmic (T) spectral density. In 
addition, they are compared with the approximate analytical, weak- 
coupling results of Eqs. i 1 1 i and 1 12t for the superohmic (solid line) 
and the ohmic (dashed line) cases. The parameters are: (i) Ohmic 
case, q = 0.05 and (ii); superohmic case, aso = 0.075. For all 
cases, T = 15 mK. 
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FIG. 7: Q-factor as a function of the dot radius a for two different 
tunneling amplitudes A and for a fixed ratio d/a = 3. The remaining 
parameters are the same as in Fig. [2] 
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FIG. 6: Q-factor as a function of temperature for two different tun- 
neling amplitudes A. The remaining parameters are the same as in 
Fig.0 



